Experimental Requirements to Determine the Neutrino Mass Hierarchy Using 

Reactor Neutrinos 
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This paper presents experimental requirements to determine the neutrino mass hierarchy using 
reactor neutrinos. The detector shall be located at a baseline around 58 km from the reactor(s) to 
measure the energy spectrum of electron antineutrinos (V e ) precisely. By applying Fourier cosine 
and sine transform to the L/E spectrum, features of the neutrino mass hierarchy can be extracted 
from the lAmjjjJ and |Amf 2 | oscillations. To determine the neutrino mass hierarchy above 90% 
probability, requirements to the baseline, the energy resolution, the energy scale uncertainty, the 
detector mass and the event statistics are studied at different values of sin 2 (2^13). 

PACS numbers: 13.15.+& 14.60.Pq, 14.60.Lm 



Neutrino physics has undergone a revolution over the 
last decade and reaches now to an era of precision mea- 
surement for the neutrino oscillation parameters. How- 
ever, 6*i3, CP-violating phase and the sign of Am| 2 (mass 
hierarchy) are still undetermined. Usually normal hier- 
archy (NH) is defined as Am§ 2 > and inverted hierar- 
chy (IH) is defined as Am|n < 0. Accelerator neutrino 
experiments such as Nova [l|, 0, Q and T2KK Q have 
the potential to determine the mass hierarchy using the 
matter effect of neutrinos at long baselines. There are 
also discussions to precisely measure the distortions of 
the reactor neutrino energy spectrum at an intermediate 
baseline (40 — 65 km) 0, [y] . A Fourier transform method 
was recently proposed to enhance and visualize the fea- 
tures of mass hierarchy in the frequency (Am 2 ) spectrum 

A new study L 8| based on Fourier transformation uti- 
lizing both the amplitude and phase information is pre- 
sented recently to enhance the features distinguishing the 
mass hierarchy at a very small sin 2 (26*13) value. This pa- 
per is complimentary to Ref. Q by taking into account 
experimental details, such as the baseline, detector re- 
sponse including energy resolution, energy scale uncer- 
tainty, and the event statistics, etc. The study is based 
on Monte Carlo simulation. 

Taking into account the detector response, the reactor 
neutrino V e L/E spectrum F(L/E) becomes 

F(L/E') = J R{E,E')F{L/E)AE 1 



F{L/E) = ct>(E)a(E)P ee (L/E), 



(1) 



where L is the baseline, E is the actual v e energy, E' is 
the observed ~D e energy taking into account the detector 
response, and R(E, E') represents the detector response 
including effects such as the energy resolution and en- 
ergy scale. The reactor neutrino flux, 4>{E), the neutrino 
inverse beta reaction cross section with detector, cr(E), 
and the neutrino oscillation probability, P ee (E), have all 
been described in Ref. 8]. Here for completeness, we 
rewrite the V e survival probability -PeeC^)!^ as: 
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cos 2 (# 12 )sin 2 (2Msin 2 (A 3 i), 
sin 2 (0 12 )sin 2 (20 13 )sin 2 (A3 2 ), (2) 



The analytical formulas for Fourier cosine and sine 
transform [8J] are: 



FCT(lu) = 
FST(uj) = 



F(t) cos(ujt)dt, 
F(t) sm{wt)dt, 



(3) 



where u> is the frequency defined as 2. 54 Am 2 ; t = L/E is 
the variable in L/E space, varying from t m i n = L/E max 
t° tmax ~ L/E rnin . In real experiments with a set of 
discrete events, the integral can be changed to the sum- 
mation over all events as in the following: 

N 



FST(u) = J2M"L/El), 

i=l 
N 

FCT{uj) = cosiuL/El), 



(4) 



where E[ is the measured energy of individual events, 
and N is the total number of events collected. 

The actual experimental measurements of the neutrino 
energy usually have two aspects of detector responses: 
energy resolution and energy scale. The response of 
the detector due to energy resolution can usually be 
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described by a Gaussian function 
where erg is the energy resolution, 
energy are usually measured by scintillators, the energy 
is typically proportional to the number of photoelectrons, 
and the error is dominated by the photoelectron statis- 
tics. Therefore the neutrino energy resolution is propor- 
tional to f / \/E V is, where E V i s = E„ — 0.8 MeV is the neu- 
trino visible energy in the detector. Previous experiments 
typically have an energy resolution of about f 0%/\/E V i S . 
Different detectors may have different forms of the energy 
scale uncertainty. For simplicity, we take two possible 
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TABLE I: Default values for neutrino oscillation parameters 
and other input parameters studied in this paper. 



cases, shift and shrinking/expanding. It is modelled as 
formula E' = (1 + a)E + b, where a and b are parameters. 

In this study, each Monte Carlo experiment generates 
a set of v e events by sampling F(L/E) spectrum with in- 
put parameters {sin 2 #13, L, oei a, b}. The total number 
of generated events determines the statistical error. De- 
fault oscillation parameters are taken from Ref. Q and 
reproduced here in Table HI together with default input 
parameters to be studied in this paper. 

In the following study, effects of input parameters to 
the mass hierarchy are studied one by one, while the rest 
remain at the default values. It is known that at the os- 
cillation maximum of Am 2 2 , corresponding to a baseline 
of about 58 km, the sensitivity to the mass hierarchy is 
maximized. Hence the default value of the baseline is 
set at 58 km. The mass hierarchy effects in the reac- 
tor neutrino energy spectrum are mainly characterized 
by Am^/IAm^l, which is only about 3%. Hence the 
required energy resolution shall be at the level of this 
number and the default value in this study is set to be 
3% / \J E V i S . Previous experiments show that a and b are 
typically at the level of 1% and 0.01 MeV respectively. 
The observed neutrino event number is proportional to 
the detector volume, exposure time and reactor(s) power. 
A very powerful reactor complex can consist of 8 reactor 
cores, each with ^3 GW thermal power. With a base- 
line of 58 km from such a reactor complex and taking 
into account the oscillation probability which is around 
sin 2 2#i2, 5 x 10 5 events, corresponding to a detector ex- 
posure of ~700 kt-year, is taken as the default for an 
experiment. 

Fig. [T] shows FCT and FST spectra from a Monte 
Carlo simulation using parameters (sin 2 (2^13), erg) = 
(0.02, 3%/^Eyis). For comparison, the analytical spec- 
tra are also shown at sin 2 (26*13) = 0.02. The impacts 
of the energy resolution and statistical errors are obvi- 
ously seen as that the amplitudes of noisy peaks and 
valleys appear to be higher in the frequency range of 
2.0 x 10~ 3 eV 2 < Am 2 < 2.8 x 10" 3 eV 2 . However, the 
main peak and valley are distinctive and can still be used 
to determine the neutrino mass hierarchy. 

We introduce parameters RL and PV [8| to quantify 
the features of FCT and FST spectra. 
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(a)FCT spectrum. 
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(b)FST spectrum. 

FIG. 1: FCT and FST spectra from simulation with parame- 
ters (sin 2 (26»i 3 ),CT B ) = (0.02,3%/v / S«i 3 ), together with the an- 
alytical spectra for sin 2 (20i3) = 0.02. Solid and long-dashed 
lines are spectra based on simulation for NH and IH cases, 
while dashed and dotted lines are analytical spectra. 



amplitude of the peak and V is that of the valley in the 
FST spectrum. 

For each set of input parameters 
{sin 2 6*i3, L, <je, a, b, N}, we simulate 500 experiments 
and calculate the probability to determine the mass 
hierarchy based on the distributions of RL and PV 
values. The procedure is concluded as the following: 

1. Given sin 2 #13 and L, we sample V neutrino events 
with energy Ei(i = 1,2, N) from energy spec- 
trum both for NH and IH cases. 

2. Ei is smeared and/or shifted to E[ based on the 
given energy resolution (<te) and energy scale un- 
certainty (a and b) parameters. 

3. FCT and FST spectra are calculated using Eq. |4j 



where RV is the amplitude of the right valley and LV 
is that of the left valley in the FCT spectrum; P is the 



4. RL and PV values are calculated based on FCT 
and FST spectra using Eq. [S] 
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5. Repeat the above steps 500 times and obtain the 
distributions of RL and PV values. 

6. Calculate the probability to determine the mass hi- 
erarchy correctly based on the distributions of RL 
and PV values. 

Fig. [2] shows the distribution of RL and PV 
values for 500 experiments with input parameters 
(sin 2 (26>i 3 ), o E ) = (0.02,2%/VE~7). Two clusters of 
points in the (RL, PV) plane corresponding to NH and 
IH cases show the probability to determine the mass hi- 
erarchy. Various input parameters have been tried and 
the distribution of RL + PV is shown in Fig. [3] Two 
clusters of points turn into two Gaussian distributions 
and the probability to determine the mass hierarchy can 
be correctly calculated. 
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FIG. 2: Distribution of RL and PV values from 
simulated experiments with parameters (sin 2 (2#i3), <je) = 
(0.02, 2%/\/E V i 3 ). Two clusters of points show the sensitivity 
to determine the mass hierarchy. 



To study the impact of the baseline, a total of 500 ex- 
periments have been simulated for each set of input pa- 
rameters (sin 2 (26*13), C-b)- Fig. 0] shows the results. The 
error bars are due to statistics since only a limited num- 
ber of experiments are simulated. The optimal baseline 
is clearly 58 km, which is chosen as the default baseline. 

Fig. [5] shows the impact of event number to the de- 
termination probability. Obviously, fewer number of 
events will induce larger statistical fluctuations, more 
noisy peaks and valleys in the FCT and FST spectra 
and hence reduce determination probability. As shown 
in Fig. [5J a total of 5 x 10 5 events will reach 90% deter- 
mination probability for sin 2 (26 i i3)=0.02 with an energy 
resolution of 3%/y/E V i S . This number of events, proba- 
bly the largest can be imagined nowadays, is chosen as 
the default in this paper. 

The requirement to the event statistics strongly de- 
pends on the value of sin 2 (2£q 3 ). Fig [6] shows the num- 
ber of neutrino events needed to determine the mass 
hierarchy at the 90% confidence level as a function of 
sin 2 (26'i3). Two cases of the energy resolution, 2% / y/ E V i s 
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FIG. 3: Distribution of RL + PV values from 1000 
experiments for parameters (sin 2 (2(9i3), oe) being 
(0.01, 2%/-s/M~), (0.01,3%/ 'VE~), (0.02, 2%/^E~) 
and (0.02, 3%/V E V i a ). Two Gaussian distributions show the 
sensitivity to determine the mass hierarchy. 
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FIG. 4: Impact of the baseline to the determination proba- 
bility for 4 sets of parameters (sin 2 (26*13), cb). 



and 3%/i/Evis, are studied. If sin 2 (2^13) happens to be 
more than 0.05, as some of the recent global fit indi- 
cated 10], the number of events can be a factor of 5 
smaller than that in the case of sin 2 (20 13 )=O.O2. 

Impact of the energy resolution to the mass hierarchy 
determination is studied for the cases of sin 2 (2^13) = 
0.02, 0.01 and 0.005 as shown in Fig. To achieve 
the mass hierarchy determination probability better than 
90% at sin 2 (2#i 3 ) = 0.02, the energy resolution shall be 
better than 3%\/E vis . This is actually a very stringent 
requirement, at least a factor of two better than that of 
the existing reactor neutrino experiments. For a typi- 
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FIG. 5: Impact of the event number to the determination 
probability in six sets of parameters (sin 2 (2#i3), ere) 



in 

g 10 3 

0) 



S3 

E 

3 



10 2 r 



resolution = 2%A/E 
resolution = 3%A/E 




0.01 



0.02 0.03 



0.04 0.05 

sin 2 29 13 



FIG. 6: Requirements to the number of events to deter- 
mine the mass hierarchy at 90% probability as a function of 
sin 2 (20i 3 ). 



which results in the frequency spectra (FCT spectra) left 
shifted. The energy scale uncertainty only introduces a 
bias to the oscillation frequency and hence Amf^ (shown 
as the main peak in the FCT spectrum). Since our 
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FIG. 7: Impact of the energy resolution to the determination 
probability for sin 2 (20i 3 ) = 0.02, 0.01 and 0.005 
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FIG. 8: Impact of sin 2 (2613) to the determination probability 
in four cases of the energy resolution. 



cal liquid scintillator experiment, substantial more light 
shall be collected to reach such a level. 

Fig. [5] shows the impact of sin 2 (26>i3) to the determi- 
nation probability in four cases of the energy resolution 
{0//E~ 1%/vfe, 2%/VE~7, 3%/VE~}. It shows 
that the energy resolution is very important to deter- 
mine the mass hierarchy, while a larger sin 2 (2# 13 ) can 
relax substantially such a requirement. 

The impact of the energy scale uncertainty is stud- 
ied by transforming the sampled neutrino energy E to 
E' = (1 + a)E + b. For two cases of a = —1% or 
b = —0.01 MeV, which correspond to the shrinking or 
left-shift of the neutrino energy spectrum, the FCT spec- 
tra are calculated and shown in Fig. [9] It shows that the 
FCT spectra, both for NH and IH cases, are left shifted. 
After shrinking the energy spectrum, the L/E spectrum 
expands and the oscillation frequency becomes smaller, 



method only depends on the relative position of peaks 
or valleys in FCT and FST spectra, the mass hierarchy 
determination is not affected by the energy scale uncer- 
tainty. 

In summary, we have studied experimental require- 
ments to determine the mass hierarchy using Fourier 
cosine and sine transform to the reactor neutrino L/E 
spectrum. The parameters RL and PV are defined to 
extract features of the Fourier sine and cosine spectra, 
and the mass hierarchy can be determined from events 
collected in experiments similar to that in the analytical 
case. The impacts of baseline, event statistics, energy 
resolution and energy scale uncertainty to the mass hi- 
erarchy determination are studied in detail. This paper 
provides a guidance to the design of the experiment to 
determine the mass hierarchy using reactor neutrinos. 
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FIG. 9: FCT spectram for an energy shift of 0.01 MeV and 
shrinking of 1%. 



